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AND  MILLIMETER  WAVT  INTEGRATED  CIRCUITS  AND  SUBSYSTEMS 


ABSTRACT 


This  interim  technical  report  presents  results  of  research  on  the  computer  aided  design  of 
monolithic  microwave  and  millimeter  wave  integrated  circuits  and  subsystems.  A  specific 
objective  is  to  extend  the  state-of-the-art  of  the  computer  aided  design  (CAD)  of  the  monolithic 
microwave  and  millimeter  wave  integrated  circuits  (MIMIC).  In  this  reporting  period,  we  have 
derived  a  new  model  for  the  high  electron  mobility  transistor  (HEMT)  based  on  a  nonlinear 
charge  control  formulation  which  takes  into  consideration  the  variation  of  the  2DEG  distance 
offset  from  the  heterointerface  as  a  function  of  bias.  Pseudomorphic  InGaAs/GaAs  HEMT 
devices  have  been  successfully  fabricated  at  UCSD.  For  a  1  jam  gate  length,  a  maximum 
transconductance  of  320  mS/mm  was  obtained.  In  cooperation  with  TRW,  devices  with  0.15  pm 
and  0.25  pm  gate  lengths  have  been  successfully  fabricated  and  tested.  New  results  on  the  design 
of  ultra-wideband  distributed  amplifiers  using  0.15  pm  pseudomorphic  InGaAs/GaAs  HEMT’s 
have  also  been  obtained.  In  addition,  two-dimensional  models  of  the  submicron  MESFET's, 
HEMT’s  and  HBT’s  are  currently  being  developed  for  the  CRAY  X-MP/48  supercomputer. 
Preliminary  results  obtained  are  also  presented  in  this  report. 
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1.  Fabrication  of  HEMT  Devices 
1.1  Introduction 

High  electron  mobility  transistors  (  HEMT  )  have  shown  great  potential  for  high  speed  and 
high  frequency  applications.  Much  progress  has  been  made  in  the  AlGaAs/GaAs  HEMTs  where 
the  larger  band  gap  highly  doped  AlGaAs  supplies  electrons  to  the  two  dimensional  electron  gas 
(2DEG)  in  the  smaller  band  gap  undoped  GaAs.  Recently,  there  has  been  increasing  interest  in 
pseudomorphic  InGaAs/GaAs  HEMTs1'2,3  because  of  the  higher  mobility  of  InGaAs  over  GaAs 
and  the  elimination  of  the  persistent  photoconductivity  problem  associated  with  the  AlGaAs. 

Pseudomorphic  InGaAs/GaAs  HEMT  devices  with  a  gate  length  of  2  microns  and 
transconductance  of  120  mS/mm  have  been  fabricated  at  UCSD.  Using  the  fabrication  facilities 
and  process  at  TRW,  a  one  micron  gate  length  InGaAs/GaAs  HEMT  with  a  maximum 
transconductance  of  320  mS/mm  was  successfully  made.  These  results  are  comparable  to  those 
found  in  a  1  micron  pseudomorphic  InGaAs/AlGaAs  HEMT.  4  These  HEMTs  have  exhibited  a 
transconductance  of  540  mS/mm  in  a  0.1  micron  device 5. 


Material 

The  material  system  for  these  HEMT  devices  is  pseudomorphic  In  15Ga  g5As/GaAs  on 
GaAs  substrate,  grown  by  molecular  beam  epitaxy  at  UCSD.  The  advantages  of  this  material 
over  the  conventional  AlGaAs/GaAs  is  the  higher  electron  mobility  and  peak  electron  velocity  in 
the  InGaAs  and  the  avoidance  of  problems  associated  with  the  deep  donor  -  vacancy  complexes 
known  as  DX  centers  in  the  AlGaAs.  These  cause  the  threshold  voltage  to  be  dependent  on 
temperature  and  illumination  and  also  lead  to  a  persistent  photoconductivity  at  low  temperatures. 
The  lattice  mis-match  between  the  GaAs  and  InGaAs  would  normally  create  dislocations  at  the 
interface  but  by  making  the  InGaAs  layer  thin,  the  mis-match  can  be  taken  up  by  elastic  strain. 
Figures  1  and  2  show  the  material  structures  of  the  HEMTs  reported  here. 

*  T.E.  Zipperian,  Lit.  Dawson,  C.G.  Osbum,  and  I.J.  Frio.  IEDM  Technical  Digest.  1983.  pp.  696-698 

*  J.J.  Rosenberg.  N.  Benlamri.  PD.  Kirchner,  J.M.  Woodall,  and  GD.  Pettit.  IEEE  Electron  Device  Letters.  EDL-6,  pp.  491-493, 
1985. 

^  T.E.  Zipperian  and  TJ.  Dnimmond,  Electronic  Letters,  vol  21.  pp.  823-824,  1985. 

4  A.  Kctterson,  et  aL  IEEE  Electron  Device  Letters.  EDL-6,  no.  12,  pp.  628-630,  1985. 

^  P.C.  Chao,  et  aL  IEEE  Electron  Device  Letters  EDL-8,  no.  10,  pp.  489-491,  1987. 
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Figure  1 


Figure  2 


Dev  ice  Fabrication 


The  HEMTs  fabricated  were  split  source  devices  with  various  gate  lengths  and  widths.  A 
new  mask  set  was  made  for  2  micron  discrete  field  effect  transistors.  It  was  determined  that  2 
micron  geometries  was  the  minimum  dimension  that  could  be  reliably  fabricated  with  the  current 
HTG  aligners  at  UCSD.  The  devices  would  serve  as  an  initial  step  in  the  development  toward 
smaller  dimension  devices  and  be  useful  as  prototype  devices  for  investigating  various  material 
systems.  The  total  gate  width  was  300  microns  and  the  source  to  drain  spacing  was  6  microns, 
allowing  a  2  micron  alignment  tolerance  in  the  mask  set.  The  overall  dimensions  are  shown  in 
Figure  3a.  In  addition,  a  few  1  micron.  2.5  micron  and  3  micron  gate  length  devices  with  source 
to  drain  spacings  of  5  micron.  7.5  micron  and  9  microns,  respectively,  were  included  as  test 
structures.  The  1  micron  gate  length  device  mask  set  given  us  by  TRW  and  shown  in  Figure  3b 
and  3c,  has  two  split  source  devices  with  source  to  drain  spacing  of  3  microns  and  gate  widths  of 
300  and  100  microns,  respectively. 

The  HEMT  fabrication  process  at  UCSD  consists  of  three  sections,  which  are  mesa 
isolation,  ohmic  contact  and  gate.  The  mesa  isolation  is  accomplished  by  a  wet  etch  consisting  of 
phosphoric  acid,  hydrogen  peroxide  and  water  in  a  1:1:25  mixture.  For  the  ohmic  contacts  a 
AuGe  alloy  with  a  thin  Ni  cap  layer  is  used.  Furnace  annealing  is  performed  at  450°C.  Much 
effort  has  been  focused  on  the  gate  process,  which  is  the  most  critical  portion  of  the  fabrication. 
The  gate  lithography  involves  fine  geometries  and  alignment  and  the  gate  recess  controls  the 
threshold  voltage.  Thermally  evaporated  aluminum  is  the  gate  metal.  An  outline  of  the  UCSD 
process  is  attached.  The  intention  was  to  establish  a  consistent  2  micron  gate  process  followed  by 
work  on  a  1  micron  gate  process  at  UCSD.  Toward  this  goal,  there  have  been  beneficial 
opportunities  to  visit  and  use  TRW  facilities  in  Redondo  Beach,  California. 

A  two  micron  photolithography  process  was  developed  for  the  new  devices.  An  HTG 
aligner  with  mid  -  UV  (300-320  nm)  wavelength  source  was  used.  It  became  apparent  that 
contact  between  the  wafer  and  the  mask  is  the  most  critical  factor  to  achieving  uniform, 
consistent  and  accurate  resolution  of  the  mask  features.  The  best  results  come  with  use  of  vacuum 
contact,  a  chuck  designed  for  the  appropriate  sample  size,  and  a  flat  sample  backside.  As  a  result. 
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the  backside  of  the  sample  needed  to  be  cleaned  and  lapped  to  remove  the  In  residue  left  after 
MBE  growth.  Difficulties  arise  with  small  piece  samples  because  good  contact  is  harder  to 
achieve.  Among  the  principle  causes  are  mask  bow  and  non-uniform  resist  thickness.  Mask  bow 
can  be  reduced  with  a  mask  holder  and  chuck  designed  for  small  pieces. 

The  alignment  tools  and  microscope  optics  on  the  HTG  are  inadequate  for  alignment  within 
1  micron.  Consistent  resolution  of  1  micron  spaces  has  been  difficult.  The  Karl  Suss  mask 
aligner  is  designed  for  use  in  the  1  micron  and  submicron  range,  and  is  more  widely  accepted  and 
commonly  used.  The  fabrication  facility  at  TRW  has  Karl  Suss  aligners  and  these  were  used  for 
the  1  micron  InGaAs/GaAs  HEMT  devices  fabricated  there. 

The  gate  recess  etch  is  the  most  critical  step  in  the  processing  of  HEMT  devices.  It  is 
necessary  for  the  removal  of  the  conduction  path  in  the  highly  doped  top  layer  and  to  adjust  the 
threshold  voltage.  As  a  result,  there  was  much  development  for  this  part  of  the  process.  A 
method  was  needed  to  determine  the  endpoint  of  the  etch.  This  was  accomplished  by  monitoring 
the  saturation  current  between  the  source  and  drain  by  making  contact  to  the  ohmic  pads  through 
the  resist.  The  sample  would  be  etched  for  a  certain  time  and  measured  until  a  pre-determined 
endpoint  saturation  current  was  reached.  Two  GaAs  MESFETs  and  a  InGaAs/GaAs  HEMT  with 
five  micron  gate  lengths  were  fabricated  as  trial  runs  for  using  this  technique.  The  results  of 
these  devices  are  seen  in  Figures  4,5.  and  6.  The  gate  recess  allows  the  pinchoff  in  the  depletion 
mode  MESFET  and  the  fabrication  of  enhancement  mode  devices.  The  depth  of  the  gate  recess 
may  only  be  a  couple  of  hundred  angstroms  so  a  very  slow  etch  is  necessary. 

As  a  result,  the  next  issued  addressed  was  to  improve  the  wet  chemical  etch  system.  An 
etch  of  phosphoric  acid,  hydrogen  peroxide  and  water  in  a  1:1:25  ratio  had  been  used.  However, 
the  etch  rate  of  30-40 A/sec  was  to  fast  for  a  shallow  gate  recess  and  resulted  in  a  lack  of  control 
and  uniformity.  The  citric  acid  etch  system  was  investigated  by  Otsubo6  and  found  to  be  an 
excellent  stow  and  linear  etch  of  GaAs.  A  mixture  of  50%  aqueous  citric  acid  with  hydrogen 
peroxide  in  a  40:1  ratio  proved  to  have  an  constant  etch  rate  of  10A  sec.  This  new  etch  provided 
more  stability  and  by  increasing  the  etch  time  improved  the  control  of  the  gate  recess  depth. 

Visits  to  TRW  this  summer  and  the  opportunity  to  work  in  a  well  equipped  facility  with  an 
established  process  was  very  helpful.  Fabrication  of  InGaAs/GaAs  HEMT’s  at  TRW  serve  as  a 
benchmark  for  improvements  of  the  process  at  UCSD  and  help  qualify  portions  of  the  process  as 
well  as  identify  problem  or  critical  areas.  The  2  micron  HEMT  devices,  #348  A  and  #462C 
shown  in  Figures  7  and  8  had  the  gate  portion  of  the  process  done  at  TRW  using  the  resist  and 
resist  process  from  UCSD.  The  results  indicated  that  the  mesa  isolation  and  the  ohmic  contact 
process  as  well  as  the  gate  resist  process  are  adequate  to  successfully  fabricate  2  micron  devices 
at  UCSD.  The  complete  fabrication  process  was  then  done  at  UCSD  resulting  in  working  devices 
on  #348C  and  #462D.  However,  as  seen  in  Figures  9  and  10,  the  performance  of  these  devices 
were  not  as  good  as  those  finished  at  TRW,  indicating  some  improvement  is  needed,  particularly 
in  the  gate  metallizations. 

Sample  #462E  and  #348D  were  fabricated  with  the  TRW  ohmic  and  gate  process.  The 
results  from  the  2  micron  devices  on  #462E  were  compared  to  devices  either  partially  or 


M.  Otsubo.  T.  Oda.  H.  Kumabe,  and  H.  Miku/  Eltctrochtmical  Soc.  134,  155  (1977). 


completely  fabricated  at  UCSD.  The  performance  of  these  devices  on  #462E  and  the  1  micron 
devices  on  #348D  serve  as  a  target  for  future  devices. 

Results  and  Discussion 

The  best  results  for  the  pseudomorphic  InGaAs/GaAs  HEMTs  are  from  the  1  micron 
devices  on  sample  #348D.  A  short  gate  device  (Figure  1 1)  (gate  width  of  100  microns)  exhibits  a 
transconductance  of  320  mS/mm.  while  the  best  value  for  the  long  gate  device  Figure  12  (gate 
width  of  300  microns)  is  225mS/mm.  The  1-V  curves  shown  in  Figures  11  and  12  show  both 
positive  and  negative  gate  voltages  with  the  peak  transconductance  occurring  with  positive  gate 
voltage.  They  also  show  good  saturation  and  pinch-off  characteristics  as  well  as  good  current 
values.  These  results  are  comparable  to  the  1  micron  InGaAs/AlGaAs  pseudomorphic  HEMT  by 

4 

Ketterson  et.  al.  which  had  a  transconductance  of  270  mS/mm.  The  potential  of  the  device  is 
shown  here  and  further  improvements  to  the  process  should  increase  its  performance.  High 
frequency  measurements  will  be  made  on  these  devices. 

For  the  2  micron  devices,  the  highest  transconductance,  183mS/mm  for  a  300  micron  wide 
gate,  was  from  sample  #462E  which  was  fabricated  with  the  TRW  process  (Figure  13).  As  a 
comparison,  #462C  (Figure  8)  and  #348A.  (Figure  7)  which  only  had  gate  recess  and 
metallization  done  at  TRW,  have  transconductance  of  135  mS/mm  and  175  mS/mm.  It  was 
found  that  the  ohmic  contact  resistance  for  #462E  is  half  of  that  for  #462C.  Devices  from  the 
same  material  but  completely  fabricated  at  UCSD,  #462D  (Figure  10)  and  #348C  (Figure  9),  have 
transconductances  of  80mS/mm  and  120mS/mm.  With  some  further  development  of  the  gate 
recess  and  metallization  the  performance  of  the  devices  fabricated  at  UCSD  will  improve  and  a  2 
micron  process  will  be  finalized.  This  gives  us  the  capability  of  initial  evaluation  of  different 
material  systems  and  with  a  few  improvements  to  the  processing  equipment,  allow  us  to  establish 
a  1  micron  process  at  UCSD. 
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1.  APPENDIX  A 
UCSD 

Mesa  Isolation 
Lithography 

Solvent  clean 

Prebake  (  105  C.  20  min  ) 

Hunt  HPR  1182  photoresist 
Spin  (5000  rpm.  20  seconds) 

Softbake  (105  C.  30  min) 

Exposure  HTG  aligner  C40mJ/cm:) 
Develop  (  1  min,  no  agitation) 

Hardbake  (1 10  C.  30  min) 

Etch  (45  sec,  H3P04:H202:H20.  ratiol  :1:25  ) 
Resist  strip  in  acetone 
Ohmic  Contact 
Lithography 

(same  as  above) 

Metal  Evaporation 

AuGe:Ni  (1000 A.  thermal  evaporation 
Liftoff,  in  acetone 
Alloy 

450*  C,  30  sec  in  formine  gas  using 
furnace  with  graphite  boat 


Gate 

Lithography 

(Same  as  above) 

Recess  etch 

Preclean  (NH»OH:H;0.ratiol:15) 

Etch  ( citric  acid  :  H202:H20,  ratio  20:1:20  ) 
Measure  source  to  drain  current 
Metal  evaporation 

Preclean  (  NHiOH:H20,  ratio  1:15  ) 

O 

A1  (1000 A  thermal  evaporation 
Liftoff  in  acetone 


UCSD  348C  InGaAs/GaAs  HEMT  with  gm  =  120mS/mm 
Figure  9 


UCSD  462D  InGaAs/GaAs  HEMT  with  gm  =  80mS/rmn 
Figure  10 


UCSD  348D  InGaAs/GaAs  1pm  HEMT  (gate  width  =  100  urn) 
gm  =  320mS/mm 
Figure  1 1 


UCSD  348D  InGaAs/GaAs  1  jam  HEMT  (gate  width  =  300  (tm) 
gm  =  225mS/mm 
Figure  12 


2.  Design  of  Distributed  Amplifiers  I'sing  0.15  pm  Pseudomorphic  InGaAs/GaAs 
HEMT’s 

2.1  Introduction: 

Distributed  amplifiers  using  0.15  pm  pseudomorphic  InGaAs/GaAs  HE.MT's  for  60  GHz 

and  possibly  94  GHz  frequency  bands  are  being  designed  at  UCSD.  The  potential  of  distributed 

1 

amplification  for  obtaining  gain  over  wide  frequency  bands  has  long  been  recognized.  ’~ 
Broadband  distributed  amplifiers  using  discrete  GaAs  MESFET's  have  been  demonstrated  at 
microwave  frequencies.^'  Better  performance  is  expected  with  the  use  of  HEMT  due  to  its 

supenor  high  frequency  characteristics  over  MESFETs.  The  principle  behind  distributed 
amplifiers  was  based  on  the  idea  to  separate  the  input  capacitance  of  the  active  devices  by  means 
of  artificial  transmission  lines,  while  adding  their  transconductances.  As  a  result,  it  is  possible  to 
achieve  extremely  wide  bandwidths  up  to  very  high  frequencies  with  effectively  very  wide  gates. 
The  latter  obviously  has  a  beneficial  impact  on  the  amplifier's  power  handling  capabilities. 

For  the  distributed  amplifier  design,  we  will  first  develop  the  D.C.  model  and  small  signal 
equivalent  circuit  for  HEMT.  These  will  be  initial  steps  toward  developing  a  systematic 
approach  to  the  circuit  design  of  distributed  amplifier.  In  the  following  sections,  the  small  signal 
equivalent  circuit  will  be  discussed,  and  be  applied  to  some  distributed  amplifier  design. 


3.  Small  Signal  Equivalent  Circuit 


There  are  several  sets  of  S-parameters  measured  on  a  0.15pm  pseudomorphic  HEMT  by 
TRW.  One  set  labeled  as  HEMT  C2  is  listed  in  Table  1 .  Its  frequency  performance  MAG/MSG. 
MUG,  Si( ,  K  vs.  frequency  are  plotted  in  Fig.  1 ,  and  fm„  can  be  projected  to  higher  than  80  GHz. 
The  current  gain  h21  vs.  frequency  is  plotted  in  Fig.  2  and  fT  is  projected  to  70  GHz. 
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The  small  signal  equivalent  circuit  of  the  active  device  is  necessary  to  predict  performance 
outside  the  measured  frequency  band.  The  circuit  derived  from  the  measured  date  includes  many 
parameters,  such  as  the  parasitics  associated  with  the  fabncation  process.  Therefore,  it  is 
important  to  derive  an  exact  equivalent  circuit.  The  parasitic  parameter  resultmg  from 
measurement  should  be  de-embedded  out  of  the  small  signal  equivalent  circuit.  This  de- 
embedded  model  is  used  in  the  monolithic  amplifier  circuit  design. 

A  small  signal  equivalent  circuit  proposed  by  Dr.  Ahlgren  is  shown  in  Fig.  3(a)  to  fit  the 
measured  data  listed  in  Table  1.  and  the  deembedded  model  is  shown  in  Fig.  3(b).  An  amplifier 
with  lumped  elements  was  designed  and  is  shown  in  Fig.  3(c),  and  its  performance  is  shown  in 
Fig.  4. 

Although  the  fit  between  the  model  and  the  measured  data  is  very  good,  a  more 

conventional  model  is  derived  based  on  the  physical  structure  of  the  device  and  from  the  point  of 
8 

view  of  statistics.  In  addition,  the  parasitic  resistance  Rg  Rs  and  Rd  can  be  obtained  from  dc 
measurements  as  detailed  by  Fukui.^  A  resulting  small  signal  equivalent  circuit  is  shown  in  Fig. 
5. (a),  and  the  de-embedded  model  is  shown  in  Fig.  5.(b). 


4.  Verification  of  the  Potential  of  Distributed  Amplifier 

The  performance  of  distributed  amplifier  is  examined  to  determine  whether  its  potential 
warrants  the  development  of  a  systematic  approach  for  design.  The  de-embedded  model  derived 
by  Ahlgren  is  shown  in  Fig.  3(b),  and  its  performance  is  shown  in  Fig.  6.  Transmission  line 
sections,  whose  characteristic  impedances  are  all  50  ohms,  were  introduced  between  gate  ends 
and  between  drain  ends  of  this  de-embedded  model  to  form  a  distributed  amplifier,  which  is 
shown  in  Fig.  7.  Its  performance  within  the  frequency  band  0-30  GHz  is  shown  in  Fig.  8. 
Comparing  the  curves  in  Fig.  6  and  Fig.  8,  it  is  clear  that  the  gain  flatness  is  improved,  and  the 
return  losses  of  input  and  output  are  much  smaller. 

An  attempt  was  made  to  improve  the  performance  of  distributed  amplifier  by  using  tapered 
transmission  lines.  ^  These  transmission  line  sections  have  different  characteristic  impedances 
from  one  another  with  linearly  or  exponentially  tapered  relations.  The  circuit  topology  is  shown 
in  Fig.  9,  and  its  performance  is  shown  in  Fig.  10.  It  is  clear  that  gain  is  8.5  ±0.5  dB,  and  return 
losses  are  below  -10  dB. 
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Similarly,  using  the  de -embedded  model  shown  in  Fig.  5(b).  we  designed  another 
distributed  amplifier  with  tapered  transmission  line,  which  is  shown  is  Fig.  11,  and  its 
performance  is  shown  in  Fig.  12. 


5.  Systematic  Design  approach  to  be  Developed 


Kaiowing  that  the  distributed  amplifier  has  potential  to  increase  the  gain-bandwidth  product, 
systematic  study  should  commence  to  achieve  a  practical  design  guideline  for  distributed 
amplifiers.  These  guidelines  will  be  established  from  examining  the  first  order  performance  of 
distributed  amplifier  designs  using  the  conventional  simplified  model  rather  than  the  de- 


embedded  model.  Then  the  small  signal  equivalent  circuit  and  de-embedded  model  should  be 
re-evaluated  by  comparison  with  the  measured  S-parameters  and  dc  measurements.^  This  de- 


embedded  model  would  be  used  to  optimize  the  circuit  design  when  running  the  Touchstone 


program. 
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3.  An  Analytical  HEMT  model  based  on  nonlinear  charge-control  formulation 

3.1  Introduction 

The  high  electron  mobility  transistor  (HEMT).  a  field-effect  transistor  that  takes  advantage 
of  the  novel  properties  of  the  two-dimensional  electron  gas  at  the  heterointerface  has  shown 
promising  performance  in  the  of  high-speed  and  high-frequency  regimes  since  us  demonstration 
in  1980.  Many  models  for  HEMT's  based  on  a  simple  linear  charge  control  model  similar  to  that 
for  conventional  MOSFET’s  with  only  the  gate  insulator  layer.  The  simulations  of  the  1-V 
characteristics  of  these  devices  have  shown  fair  agreement  with  experimental  data  within  certain 
operation  regions  or  in  some  respects  [l]-[4].  Recently.  Hughes  et  al..  [5]  has  pointed  out  the 
importance  of  nonlinear  effects  on  the  charge  control  model  and  considered  the  nonlinear 
variation  of  Ad  (the  distance  of  2DEG  from  heterointerface')  utilizing  a  numerical  method.  These 
results  led  to  a  better  fit  of  the  calculated  and  measured  transconductance  characteristics. 
However,  the  physical  mechanism  cannot  be  easily  recognized  in  the  derivation  of  this  model. 
Furthermore,  the  variation  of  charge  density  after  current  saturation  was  not  considered. 

An  analytic  I-V  model  of  the  HEMT,  which  can  simulate  the  dc  wuth  the  whole  operation 
range  of  the  dc  characteristics  through  the  whole  range  of  operations  is  derived  in  this  study. 
This  model  will  be  applied  to  conventional  AlGaAs/GaAs  HEMT’s  so  that  comparisons  can  be 
made  to  the  models  previously  mentioned.  In  addition,  the  AlGaAs/GaAs  HEMT  material 
system  has  been  much  studied  so  information  is  available  on  many  of  the  material  parameters 
needed  for  the  device  analysis.  This  allows  the  models  to  focus  on  representing  the  basic 
characteristics  and  principles  behind  the  operation  of  HEMT’s.  The  models  can  later  be  readily 
adapted  to  HEMT’s  of  new  material  systems,  such  as  the  pseudomorphic  InGaAs/GaAs  HEMT 
being  developed  at  UCSD,  as  soon  as  material  parameters  become  known. 

3.2  Nonlinear  Charge  Control 

The  first  proposed  conventional  charge  control  model  for  HEMT  [1]  derived  from  Poisson's 
equation  neglected  the  Fermi-potential  term. 

^  =  T~T  ~  VT  -  VF]  (1) 

di+d* 

V  *  ^  qNc"  H- 

vt  =  4>b  - - ^ ^  i 

q  2e, 


where  n,  is  the  2DEG  sheet  earner  concentration.  E|  and  d|  are  the  permittivity  and  thickness  of 
earner  supplying  layer,  respecuvely.  q  is  the  electronic  charge,  VG  is  the  gate-source  biased 
voltage,  4>b  is  the  barrier  height  of  the  Schottky  -gate,  AEC  is  the  conduction  band  discontinuity, 
N'di  is  the  doping  concentration  in  the  top  laser,  d«  is  the  spacer  layer  thickness,  and  VF  is  the 
Fermi-potential  with  reference  to  the  conduction  band  edge  of  the  2DEG  channel. 
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Drummond  et  al.  [2]  extended  this  model  to  include  the  variation  of  Fermi  potential  of  the 
two-dimensional  electron  gas  by  a  linear  approximation  with  respect  to  sheet  carrier 


concentration. 


Vp  =  Vp0  +ans 


which  yielded 


qn,=  — 7— -  [VG-VFo-VT] 
di+d^+Ad 


where 


This  model  was  further  modified  [5]  to  take  into  account  the  nonlinear  dependence  of  the 
position  of  2DEG  channel  Ad  from  heterointerface  on  gate  bias.  The  2DEG  offset  is  calculated 
from  the  weighted  average  of  all  energy  levels. 


f  V  n4n  d„ 

n 


n  =  0.  1.  2.  3, 


where 


n,  = 


n  =  0,  1,  2,  3. 


Clearly,  in  the  above  derivation  the  computations  are  very  complicated. 

According  to  the  triangular  potential  well  approximation  the  relation  between  n,  and  the 
Fermi-level  can  be  established,  only  if  the  lower  and  the  excited  subband  are  being  considered. 


=  —  ln< - -  e  +  e  u  + 


;/)  i  q>*  H  j/i 
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As  suggested  by  eq.  (7),  we  can  approximate  EF  by 


Ef  =  a(bnj 


with  a  high  accuracy.  Substituting  eq.  (7)  into  eq.  (1),  further  manipulations  yield  the  2DEG 
sheet  carrier  concentration 


where 


(d, -Kle )[1  +  £-ns  3  ] 

q 


AEC  qND,  , 

VG  ~  dl 

q  2e, 


P=-ii^-b3 

V  di+d. 
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Eq.  (9)  indicates  that  Ad  is  proportional  to  ns  ?  .  This  corresponds  to  the  results  calculated  from  a 
variational  method  by  Stem  and  Das  Sarma  [6) . 

3.3  Drain  Current  Characteristics 


When  a  doped  layer  with  wider  band-gap  is  grown  on  top  of  an  undoped  layer  with 
narrower  band-gap,  a  two-dimensional  electron  gas  will  form  at  the  heterointerface  due  to  the 
charge  transfer  driven  by  the  conduction  band  discontinuity.  In  normal  operation  mode,  the 
depletion  of  the  doped  layer  is  attributed  to  the  charge  transfer  and  Schottky  gate  depletion.  At 
thermal  equilibrium  the  undepleted  channel  width  of  parasitic  MESFET,  h,  as  shown  in  Fig.  1,  is 

-  i 


h  =  d  i 


NDi 


qNDi 


(4>b  ~  _  vg  ) 


(11) 


where  Vn  is  the  potential  difference  between  the  conduction  band  edge  and  the  Fenni-energy  in 
the  neutral  substrate. 


The  threshold  voltage  of  the  parasitic  MESFET's  conduction  where  h  =  q>  may  be  easily 
obtained,  and  can  be  written  as 


vc  =  <5>b  ~  Vn 


qND, 

2e, 


d,  - 


nso 
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(12) 


If  the  gate  applied  voltage  VG  >  Vc.  then  the  doped  layer  is  not  fully  depleted  and  the 
conduction  of  parasitic  MESFET  should  be  taken  into  account. 

Experimental  measurement  has  shown  the  velocity  saturation  phenomenon  existing  in  the 
high  electric  field  regime  for  III  -  V  compound  material.  For  simplicity  in  mathematics 
manipulation,  the  field-dependent  velocity  model  is  approximated  by  the  following  model  for 
both  the  parasitic  MESFET  Channel  and  the  two-dimensional  electron  gas  channel. 


•  l 


v(E)  = 


pE 


(1  +  i7> 


E  <EC 


(13) 


=  V), 


E  >  Er 


(14) 


where  p  is  the  low-field  mobility,  and  Ec  is  the  critical  field  at  the  onset  of  velocity  saturation. 


3.4  Normal  HEMT  Operation  Mode 


1.  Linear  Region  (VD  <  Vsal) 

When  the  applied  drain  voltage  is  not  high  enough  to  accelerate  carriers  up  to  saturation 
velocity,  the  drain  cument  can  be  evaluated. 


rv'.vv< 
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where  w  is  the  gate  width,  a  is  the  nonlinear  charge  control  coefficient  and  V  is  the  channel 
potential. 

Integrating  eq.  (15)  ID  can  be  written  in  terms  of  y0  and  yD: 
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w  here  y0  and  yD  are  defined  as 


(dj+4)  [L+- 


>'o  -  [VG  -  VT] 


>D  =  [VG  -  VT  -  VD]  3  (18) 

2.  Saturation  Region  (VD  >  Vsat) 

In  saturation,  we  make  use  of  eq.  (9)  and  (14).  The  2DEG  channel  current  in  saturation 
region  leads  to  the  equation 

w£,H;[Vg-Vt— Vsa,]Ec: 

ID  = - - - - — -  (19) 

(di+d*))  1  +  -^-a(VG  -  VT  -  VS4t)~  3  M  + 

q  E*. 

V-  J  -  “  „ 

To  ensure  continuous  current  transition  from  linear  region  to  saturation  region,  we  make  use 
of  the  current  continuity  between  these  regions  and  obtain  the  following, 


-6(ys+£ot)<  (£oc)*-"  0-?-yS)]  +  (£5L)6  ln[ 


ys+^-a 


n  q 


q  v 


=  >s  Ec;L(l-K)  +  VG -VT -y53 


Channel  length  modulation  factor  K  and  vs  are  defined  as 
L  -  Lc 

tr  v 


y5=  vG-vT-vs 
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We  might  determine  K  and  Vsa(.  by  solving  a  set  of  two-dimensional  Poisson  equations  by 
Laplace's  method  and  keep  the  lowest  space  harmonic  [7],  However,  a  simpler  one-dimensional 
approximation  to  Poisson  equation  is  adopted  in  the  velocity  saturation  region,  which  assumes 

w  here  the  variation  of  carrier  concentration  in  the  direction  perpendicular  to  the  current  flow  is 
assumed  to  be  negligible. 

Solving  eq.  (231  subjected  to  boundary  conditions.  E(Lc)  =  El:  and  V(Lc)=Vsal,  the 
voltage  drop  across  the  saturated  channel  region  can  be  written  in  terms  of  K  and  Vsat. 


V,al  =  VD  - 


(di-t-deH  l+-‘--a(VQ— V-[~VsaI)  1+— — 


-Er,LK 


Note  that  Eq.  (20)  and  Eq.  (24)  are  functions  of  two  dependent  variables,  ys  and  K.  Thus  Vsal  and 
K  can  be  obtained  by  solving  this  set  of  coupled  equations.  Even  though  the  equations  are  of 
high-order  functions  of  ys,  the  range  in  which  the  roots  ys  and  K  are  located  is 

0  <  K  <  1  (25) 

-  -3-a  <  ys  <  y0  (26) 

P 

This  will  restricts  the  number  of  reasonable  roots  to  one. 

Once  the  gate  voltage  becomes  larger  than  Vc ,  a  conduction  channel  forms  in  the  parasitic 
MESFET,  i.e.,  the  carrier  supplying  layer  is  not  fully  depleted,  causing  to  current  contribution 
from  2DEG  channel  as  well  as  a  MESFET  channel  underneath  the  Schottky  gate. 

3.5  Parasitic  MESFET  Conduction  Mode  VG  >  Vc 

Consider  the  case  where  the  MESFET  channel  is  pinched  off  at  x  =  LM  by  the  gate  bias. 
The  channel  potential  at  this  point  reads 

VM=Vp-(0B-Vn)  +  Vc  (27) 


where  the  pinch-off  voltage  is  given  by 


q.ND1  d|  - 


Within  the  conduction  channel,  the  2DEG  will  not  be  perturbed  by  gate  and  drain  voltages, 
where  the  channel  potential  is  lower  than  VM,  and  will  i.iy  at  thermal  equilibrium  values,  n,c, 
which  can  be  estimated  self-consistently  by 
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1.  Linear  NtESFET  Channel  i  VD  <  \\, » 

The  current  flowing  through  the  NtESFET  channel  can  be  derived  as  in  the  conventional 
NtESFET  model,  which  is  given  bv 


*qFlN'oi  d  _ 
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(30) 


Since  vD  <  vM,  the  whole  2DEG  channel  behaves  like  a  linear  resistor,  the  current  flows 
through  the  2DEG  channel  reads 


I.  = 


w£  1  M-2  1  °B  _N  r.  ~  ^  p  t]  VD 

(d,+d«)  L  +  — — 

q  Ed 


(31) 


The  total  drain  current  ID  is 
Id  =  h*  I: 


(32) 


2.  Pinched-off  NtESFET  Channel  (VM  <  vD  <vitl) 

The  2DEG  current  due  to  equilibrium  two-dimensional  electron  gas  can  be  calculated  as  in 
the  linear  region. 


<*,U:  10b-Vb-Vp-Vt1Vm 

l.d.-KUl-*-  C 3]  [L-— ^-] 
q  Ec; 


(33) 


While  the  current  flows  through  the  NtESFET  channel  it  is  evaluated  with  the  boundary 
condition,  ViLm  >  =  vM .  and  can  be  expressed  as 
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Sunilar  to  the  formulation  in  the  normal  operation  mode,  the  2DEG  current  can  be 
explicitly  expressed  as 


-3wh,£,< - - - - - - —  *  V  ( -E-ai^lvp-v^l  -t-  (  )6ln[ 
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3.  Saturated  MESFET  Channel  (V.„  <+vD> 

In  this  operation  region,  formulation  of  the  current  flowing  through  the  MESFET  channel  is 
identical  to  eqn.  (34).  And  the  2DEG  current  is  very  similar  to  that  in  the  pinched-off  region  and 
can  be  derived  as 
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Similarly,  the  two  dependent  variables,  y,  and  K.  can  be  solved  from  the  following  two 
equations,  derived  from  current  continuity  requirement,  and  from  a  Poisson  equation  in  the 
saturation  region,  separately. 
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3.6  Discussions 

The  performance  of  the  model  is  evaluated  with  the  data  from  from  Drummond  et  al.[2],  we 
have  the  model  to  evaluate  its  performance.  Fig.  2  shows  the  drain  I-V  characteristics  for  the  1 
pm  device.  The  used  parameters  are  listed  in  Table  I.  This  model  gives  a  better  agreement  with 
the  experimental  I-V  curves  than  does  the  two-piece  model.  In  this  simulation  run.  the  parasitic 
source  and  drain  resistance  are  not  included  so  parasitic  resistances  are  not  taken  into  account. 
However,  this  kind  of  gain  reduction  can  be  covered  by  adding  two  parasitic  resistances  R,  and  R„ 
to  the  equivalent  circuit  model  externally. 


In  the  deeper  gate  reverse  bias  region,  the  modeled  results  have  a  better  fit  to  the 
expenmental  data  than  to  existing  models  due  to  the  more  accurate  nonlinear  charge  control 
model  operating  in  the  saturation  region.  There  the  earner  concentration  in  2DEG  channel 
shows  high  sensitivity  to  the  vanation  of  gate  potential,  because  the  scarce  earners  are  repelled 
further  from  the  heterointerface  by  the  gate  field.  This  effect  leads  to  the  reduction  in 
transconductance. 

The  next  step  to  improving  the  accuracy  of  this  model  is  further  modification  of  the 
parasitic  MESFET  pomon.  The  existing  parasitic  MESFET  model  can  not  describe  its  external 
behavior  adequate!)  from  the  view  point  of  gate  transconductance.  The  remedy  being  planned 
will  include  some  appropnate  terms  in  the  formulation  to  take  care  of  real  space  transfer  effects. 
Those  current  components  will  be  weighted  to  guarantee  the  smooth  transition  from  2DEG 
current  dominated  region  to  the  transconductance  compression  region  around  the  onset  of 
parasitic  MESFET  conduction. 

3.7  Conclusions 

We  have  derived  I-V  characteristics  model  for  the  HEMT's  based  on  a  nonlinear  charge 
control  formulation  which  takes  into  consideration  the  variation  of  the  2DEG  distance  offset  from 
the  heterointerface  as  a  function  of  bias.  The  model  has  been  kept  in  a  simple  and  analytical 
form,  making  it  a  possible  module  to  be  implemented  in  the  circuit  simulation  program  without 
adding  too  heavy  of  a  computation  load.  Also,  the  modeled  results  are  in  fair  agreement  with 
experimental  data,  which  demonstrates  that  the  model  could  be  applied  to  the  optimization  of 
device  performance. 


Schematic  diagram  illustrating  current  in  two-dimensional 
electron  gas  channel  and  parasitic  MESFET  channel 
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Comparison  of  this  model  with  Drummond's  two-piece 
velocity  model  [2] 
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4.  SDA:  A  Highly  Efficient  Two-Dimensional  Finite-Element  Simulator  for 

Semiconductor  Device  Analysis 

4.1  Introduction 

The  further  reduction  in  geometrical  dimensions  of  today's  semiconductor  devices  and  the 
advent  of  novel  devices  employing  new  materials,  has  highlighted  the  need  for  simulation, 
especially  numerical  simulation.  This  has  become  a  powerful  tool  for  providing  insight  into  the 
functioning  of  devices  by  examining  the  distribution  of  various  physical  quantities  in  the  mierior 
of  a  device  as  well  as  external  device  behavior.  In  addition,  the  predictive  capabilities  of  detailed 
numerical  device  simulations  can  be  used  to  enhance  device  characterization  techniques  and  to 
form  a  basis  for  the  technological  optimization  of  process,  device  and  circuit  designs. 

A  new  approach  to  nonlinear  finite  element  numerical  analysis  for  semiconductor  devices, 
developed  by  the  VLSI  Research  Laboratory  at  UCSD,  has  been  implemented  in  a  general 
purpose  two-dimensional  finite  element  program  for  Semiconductor  Device  Analysis  (SDA).  In 
addition  to  the  conventional  finite-element  algorithm  constructed  to  carry  out  all  types  of  linear 
and  nonlinear  finite -element  analysis  procedures,  the  following  new  features  are  incorporated  to 
enhance  computing  capability,  accuracy,  and  graphical  interaction.  In  SDA.  mesh  generation  for 
an  arbitrary  semiconductor  device  and  the  choice  of  an  optimal  solving  algorithm  can  be  done 
automatically  by  SDA  or  interactively  through  a  dialogue  with  the  user.  By  using  high  order  error 
estimation  and  predictive  techniques,  a  strategy  for  fully  adaptive  spatial  grid  generation  and 
optimal  choice  of  the  time  steps  developed  to  find  a  balance  between  the  reduction  of 
discretization  errors  and  the  additional  amount  of  computational  efforts  needed.  The  high 
computation  speed  of  SDA  is  achieved  by  using  simultaneously  a  modified  Scharfetter-Gummel 
method  to  solve  the  two-dimensional  Poisson's  equation  and  the  two-dimensional  continuity 
equation  on  an  optimal  self-adaptive  triangular  mesh.  New  physical  models  have  been  derived  for 
the  analysis  of  and  submicron  devices  and  new  IH-V  compound  devices,  such  as  the 
pseudomorphic  InGaAs  HEMT.  These  will  be  useful  for  investigating  the  device  physics  as  well 
as  suggesting  new  approaches  for  improvement. 

This  program  can  simulate  GaAs  MESFET  and  AJGaAs/GaAs  heterojunction  bipolar 
transistors  as  well  as  conventional  silicon  MOSFET  and  bipolar  transistors.  These  devices  can 
have  arbitrary  geometry  and  impurity  distribution  with  a  wide  temperature  range. 

The  program  also  considers  incomplete  ionization  and  variable  band  structure  effects  and 
can  employ  either  Boltzmann  or  Fermi  statistics.  The  graphics  capability  makes  the  graphic 
interpretation  of  simulation  results  very  clear  and  attractive. 
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4.2  Basic  Equations 

The  two-dimensional  numerical  model,  SDA-1A.  developed  at  UCSD  for  semiconductor 
device  analysis  first  examines  Poisson's  equation 

T  u'V^i  =  n  -  p  ♦  N,  (1) 

w  here  n.  p  are  electron,  and  hole  concentration,  respectively  v  is  potential  and  N.  is  net  impurity 
concentration,  w  hich  is  defined  as 


N,  =  Nf,-N, 

Nf,  and  N,  the  ionized  donor  and  acceptor  concentration,  respectively,  are  given  by. 
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where  E0  is  the  donor  energy  level  and  E*  is  the  acceptor  energy  level.  The  electron  and  hole 
concentration  can  be  represented  as 

n  =  «vc‘-*  (5) 

p  =  n„  e*r ~  *  v'l,v'  (6) 


where  <pa  and  op  are  the  electron  and  hole  quasi-Fermi  potentials  respectively.  v„  and  vp  are 
electron  and  hole  band  parameters  respectively.  n„  is  intrinsic  carrier  density  for  reference 
material  and  thermal  potential  is  given  by 


v.  = 


kBT 


The  equations  for  the  electron  and  hole  band  parameters  are  as  follows 
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The  potential  difference  between  the  fermi  levels  and  the  conduction  or  valence  band  are  defui&J 
as 


■  '7.  = 
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Epn-E.  2  i  N„ 


k«T 


q(  y-4i  »-i  Chi -Chi,  )-E, - fn  -p-  e~*w **T 


(9) 


where  x  is  electron  affinity.  E,  is  bandgap.  N..  N.  are  effective  density-  of  states  for  the  conductor 
band  and  valence  band  respectively.  The  subscript  r  represents  the  reference  material,  and  F,  .  is 
Fermi  integral. 


The  second  major  equation  for  device  analysis  is  the  continuity  equation. 


r-  ■"  '1n  D 

V  j„  =  —  *  R 

(11) 

<n 

c  r = -  ft-  -  r 

(12) 

p  di 

where  R  is  the  recombination  rate  which  includes  SRH  bulk  and  surface  recombination.  Auger 
recombination  and  impact  ionization.  The  recombination  rate  in  turn  is  described  by 

R  =  - (n^~n  - - *  (np~n,:>  (C  jH-C„n)  -  —  (a„  |  JJ  -  ap  |Jp  | )  (13) 

r,„(n*u,  H-r„ip*-n, »  4 


and  J„.  Jp  are  the  electron  and  hole  current  densities  respectively  are  defined  as 


(14) 

jp  =  -qpMp^p 

(15) 

where  and^p  are  electron  and  hole  mobilities.  For  GaAs,  the  following  empirical  electron 
mobility  model,  which,  is  a  function  of  impurity  concentration,  temperature  and  field,  is  used 


M(Nt.  T)  =  < 


/r,  (Nt-106N^) 

Nt  +  10‘N„f 


[1  -*-NT/N„f]a  j 


T 

300 


(16) 


M(Nt.  T,  E)  = 


1  ♦  ( I E  |  /Ec  )* 


(17) 


v-  =  v.«(300Ki-/3T  (18) 

where  E  is  field  and  Vltl  is  saturation  speed. 

It  should  be  emphasized  that  the  model  described  above  has  taken  into  account  incomplete 
ionization  of  impurities  (  No  *  ND.  N;  * NA  in  (2).  (3),  (4)  ).  and  the  different  conditions  which 
require  either  Fermi  or  Boltzmann  statistics.  In  equations  (7)  and  (8),  with  Fermi  statistics 
Fi/2(bc) Fw,(tj,) and  for  Boltzman  statistics  F,/2(tjc)  =  c1c,  Fu:(Tjj  =  en'.  Finally,  for 
heterojunctions  we  find  that  v„#o,vp*oin  (5),  (6). 
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4J  Discretization  Method  Using  Scharfetter-Gummel  Scheme 

The  equations  (1),  (11),  and  (12)  as  described  above  can  be  expressed  in  a  general  form  like 

T.aVu)  =  f  (19) 


where  <t>„.  <t>p). 


f  =  (n-pi-N,.  4^— R.  - ).  a= 
F  d\  dt 


^  (£.  e' 


After  manipulation  using  a  finite  discretization  scheme,  the  discretized  equation  set  becomes 

F,  =  JaVuO  ilil  +Jf<5,dil  =  0  (20) 


This  is  a  set  of  non-linear  equations  that  can  be  solved  by  Newton  method.  The  equation  set  in 
Newton  iteration  form  is 


IlK,,  +<3,,)  Vu,  =  -  £K,,u,+f, 

'  1 1 


(21) 


where 


K,,  = 


=}-^-K»,f,dn 


chi 

f,  =J  w.dfl 

A  seven  point  Gaussian  quadrature  method  is  used  in  eqs.  (22)-(24),  which  reads  as 

|FdU  =  £  w.F, 


(22) 

(23) 

(24) 


(25) 


where  w,  is  weighting  factor  at  the  i*  point  and  F,  is  the  value  of  function  F  at  the  i*  point.  The 
finite  element  mesh  is  comprised  of  triangular  elements  and  so  the  1st.  2nd  and  3rd  points  are  the 
vertices  of  the  triangle,  the  4th.  5th  and  6th  points  are  the  mid-points  of  the  edges  and  the  7th 
point  is  the  center  of  the  triangle. 

Calculation  of  the  values  of  n.  p.  er.  e  *’  at  the  4th,  5th.  6th,  and  7th  points  is  very  important 
and  critical  because  it  has  a  significant  effect  on  the  accuracy  of  discretization  and  the 
convergence  speed.  First,  we  apply  the  Scharfetter-Gummel  (S-G)  formula  into  Poisson's 
equation  as  well  as  into  the  continuity  equations  using  finite  element  methods,  so  that  the 
computational  effort  required  for  convergence  is  reduced  and  the  accuracy  of  solutions  is 
improved  significantly. 
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The  S-G  formulas  assumed  in  SDA-IA  are  as  follows: 


n(xi  =  (  i  -  1  )n,  -  (26) 

p(si  =  ( 1  -  gix.-^li/,,  i  )p,  -  gtx.-.W,  T 
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gix.vi=  1  -  exp  y -  i  I-expiy 


Bi  X)  =  x  /  [expix  i  -  1  ] 


±V,,  =  v-w,. 


K  =x-z 


(27) 


To  illustrate  the  effectiveness  of  this  new  scheme  described  above,  which  will  be  called  method 
II,  an  abrupt  n-P  junction  is  simulated,  and  results  are  shown  as  Table  1  -  Table  4.  Method  I  is  the 
classical  finite  element  method,  using  low  order  polynomials  to  approximate  n,  p.  e*.  ande"*'. 

Table  1.  Number  of  iterations  required  for  convergence  (NUM)  versus  impurity  density  (Nd) 
with  method  II  (the  number  of  mesh  points  Np  =  24,  Na  =  1  x  to15  cm"5,  zero  bias) 


Nd  (cm-3) 

5  x  to15 

5  x  to16 

1 

5  x  to17 

1 

5  x  to1' 

5  x  to” 

1 

1  x  to35 

NP 

24 

24 

24 

24 

24 

24 

NUM 

6 

6 

7 

;  8 

8 

7 
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Table  2.  NT.N1  and  the  fewest  number  of  mesh  points  required  for  convergence  versus  Nd 
with  method  I  (classical  method)  (  Sa  =  1  x  to"  cm1,  zero  bias) 


Njicrrf '  i 

5  x  to15 

5  x  io1* 

5  x  to1’ 

5  x  io” 

5  x  to"' 

1  x  i(A 

N. 

Pr**n 

24 

37 

149 

194 

352 

446 

NT’S! 

1 1 

36 

17 

70 

134 

260 

Table  3.  A  comparison  of  discretization  errors  of  Poisson's  equation  with  method  1  and  with 
method  II.  as  a  function  of  minimum  mesh  size  aX^uh  (Na  =  1  x  to15  cm"3. 
Nd  =  5  x  to15  cm  3,  zero  bias) 


0.15 

0.075 

0.0375 

ERR 

I 

46.07 

21.30 

9.30 

% 

,  u 

13.68 

8.23 

1.58 

Table  4.  A  comparison  of  discretization  errors  of  continuity  equations  with  method  I  and  that 
with  method  II,  as  a  function  of  AX^,,  (Na  =  1  x  to”  cm*’.  Nd  =  5  x  to”  cm"3.  Vw  =  0.6V; 


!  0.15 

0.075  ; 

0.0375 

0.15  0.075 

0.0375 

ERR 

!  I  n 

34.19 

. 

19.86 

10.58  !  p  24.19 

12.21  ! 

(%) 

n 

20.20 

5.53 

2.01  i  12.01 

5.76 

From  Tables  1 

and  2,  it  is 

shown  that 

method 

0  is  sufficiently  accurate 

so  that  the 

convergence  of  the  problem  is  almost  independent  of  the  impurity  density  and  the  grid  size  of 
mesh.  This  is  in  contrast  to  method  I  where  the  convergence  is  very  sensitive  to  the  impurity 
density  and  fineness  and  layout  of  the  mesh.  This  is  due  to  large  integration  errors,  especially  in 
highly  doped  regions.  When  nd  =  to®  cm"’  and  N,  =  5xio15,  the  new  method  only  needs  7 
iterations  on  a  rather  rough  mesh  (the  number  of  points  in  the  mesh  <NP  =  24).  while  the  classical 
method  needs  246  iterations  on  a  rather  fine  mesh.  (Np  =  446).  This  shows  the  convergence  rate 
increases  by  a  factor  of  37.  The  saved  time  is  greater  because  the  time  per  iteration  on  rough 
mesh  is  far  less  than  on  fine  mesh. 
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From  Tables  3  and  4.  u  is  shown  that  with  method  II.  we  can  obiain  more  accurate  solutions 
than  with  method  I  on  the  same  mesh  particularly  tor  Poisson  s  equation.  In  addition,  the  results 
indicate  that  applying  method  II  to  the  continuity  equations  will  be  able  to  take  the  same 
advantage  of  the  S-G  scheme  as  those  well  known  in  the  finite  difference  counterpart.  When 
using  the  mesh  whose  minimum  mesh  si/e  ax^.,.  0.0375  q  the  accuracy  of  solution  ol  v.  n.  p  tor 
method  II  is  increased  by  a  factor  of  5.8.  5.26.  3.20.  respectively,  over  method  I.  The  smaller 
iv,,.  becomes  the  higher  the  improvement  of  accuracy. 


4.4  Non-Linear  Iteration  Algorithm 

It  is  well  known  that  decoupled  methods  have  the  advantage  of  higher  speed  and  less 
storage  requirements  compared  with  coupled  methods  in  the  low  and  medium  current  regime. 
However,  its  speed  is  much  low  er  than  that  of  coupled  method  in  the  high  current  regime.  This  is 
due  to  a  considerable  increase  in  the  number  of  relaxation  iterations  between  equations  needed 
for  convergence  w  hen  there  is  strong  coupling  among  the  equations.  Therefore,  we  proposed  the 
following  schemes  to  overcoming  the  drawbacks  of  decoupled  methods. 

( 1 )  Adaptive  control  of  inner  iteration  number  in  solving  a  single  equation. 

The  inner  iteration  is  New  ton  iterauon  and  the  out  iteration  is  S-G  iteration.  Because  the 
error  function  per  outer  iteration  usually  decreases  only  by  a  factor  of  03  -  0.5.  it  is  not 
necessary  to  have  the  inner  iteration  reach  convergence  during  the  initial  phase  of  outer 
iteration  befoic  the  start  of  next  outer  iteration.  Therefore,  the  following  strategy  is 
adopted: 

If  ||F J  >  C  ,  the  inner  iteration  proceeds  until  its  error  function  decreases  by  an  order  of  two. 
Otherwise  the  iteration  keeps  going  until  the  convergence  is  reached. 

The  computation  shows  that  the  desired  computauon  quantity  is  decreased  considerably  by 
the  method. 

(2)  Variable  order  prediction  technique  and  adaptive  control  of  the  step  for  bias  voltage. 

A  variable  order  and  step  method,  which  is  similar  to  the  auto-integration  algorithm,  but  in 
terms  of  voltage  step  instead  of  time  step,  is  presented.  It  is  shown  that  the  number  of  out 
relaxation  iterations  is  decreased  significantly  by  the  method. 

(3)  The  method  to  speed  up  convergence  by  using  SOR  and  extrapolation  in  outer  iteration  is 
assumed. 
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4.5  The  Self-Adaptive  Mesh  Refinement  Technique 

The  initial  mesh  for  finite  element  discretization  is  generated  by  the  simulator  according  10 
device  structures  provided  by  the  users,  and  the  self-adaptive  mesh  refinement  procedure  will  be 
invoked  to  optimize  the  mesh  whenever  the  speed-up  of  convergence  rate  is  possible.  The  seif- 
adaptive  mesh  requirement  technique  employed  in  SDA  has  eliminated  the  tedious  mesh 
specification  and  has  enhanced  the  program's  efficiency  substantially. 

4.6  Some  examples  of  Device  Analysis 

For  purposes  of  illustration,  the  SDA  program  has  been  used  to  simulate  GaAs  MESFET 
device.  The  device  dimensions  and  refined  grid  are  shown  in  Fig.  1.  Contour  maps  of  potential 
distribution  and  three  dimensional  electron  density  distribution  at  several  biases  under  thermal 
equilibrium  condition  are  plotted  in  Fig.  2  and  Fig  3.  respectively.  Fig.  4  illustrates  the  external 
I-V  characteristics  obtained  in  this  sample  run.  which  demonstrated  that  this  SDA  simulator  can 
effectively  model  semiconductor  devices  of  arbitrary  structures  within  a  wide  range  of  bias 
conditions,  and  can  satisfy  any  user  requirement. 

4.7  Further  Work 

Revisions  to  the  SDA-1A  program  are  underway  to  include  new  function  and  models,  such 
as  hot  electron  models  which  are  pertinent  to  devices  with  small  size,  and  quantum  well  and  bulk 
electron  effects  which  are  needed  for  HEMT  simulation. 
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